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We investigate the zero-temperature phases of bosons in a one- dimensional optical lattice with an 
explicit tunnel coupling to a Bose condensed particle reservoir. Renormalization group analysis of 
this system is shown to reveal three phases: one in which the linear system is fully phase-locked to 
the reservoir; one in which Josephson vortices between the one-dimensional system and the particle 
reservoir deconfine due to quantum fluctuations, leading to a decoupled state in which the one- 
dimensional system is metallic; and one in which the one-dimensional system is in a Mott insulating 
state. 



Introduction: Ultracold atoms trapped in optical lat- 
tices is an active area of research in many-body physics 
[1-6] . At zero temperature and without disorder, bosons 
in a lattice with integer filling exhibit two distinct phases 
- a superfluid phase, where the phase of the wavefunc- 
tions are sharply defined on sites, and a Mott insulating 
phase, where the site occupation numbers are sharply 
defined [7-11]. This system allows novel many-body 
phenomena associated with unusual correlated states and 
quantum phase transitions [12-19]. 

In this work we will study a system of bosons in a 
one-dimensional (ID) optical lattice that is tunnel cou- 
pled to a three-dimensional (3D) Bose-condensed particle 
reservoir. Our goal is to understand how the reservoir 
impacts the states of the bosons in the optical lattice. 
We find that the system supports three phases: a 3D 
superfluid phase in which the ID system becomes phase- 
locked with the reservoir; a decoupled phase (in which, 
under circumstances described below, the ID system be- 
haves like a metal); and a Mott insulating phase. These 
states are separated by deconfinement transitions [20], 
either of vortices or of tunneling events, as we describe 
below. The states may be distinquished both by their 
low-energy excitation spectra, and by their conduction 
properties. 

Model: Our analysis begins with a Hamiltonian for 
lattice bosons (in number-phase representation) [7] 
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where (ni T , 4>i T ) and (Nn T , $R r ) are the canonically con- 
jugate occupation number fluctuations and phases of the 
bosons in the chain and in the 3D reservoir, respectively, 



t and U > (tR and Ur > 0) are the nearest neighbor in- 
terchain hopping and on-site repulsion terms for the ID 
(3D) system, and Ji is a tunneling amplitude between 
a site on the chain [i] and a site in the reservoir [R(i)]. 
While the particle numbers in the real optical lattices are 
typically small, we nevertheless allow the fluctuations to 
vary from —00 to +00. This simplifies the calculation, 
and should not affect the allowed phases because the sym- 
metries of the resulting Hamiltonian are unchanged [21]. 
Our geometry could be realized as a linear array of sites 
using red detuned light focused near the edge of a BEC 
cloud [12,22-24]. Note that the tight-binding form for the 
reservoir is adopted purely as a matter of convenience. 

Following standard procedure [25], we construct a 
path-integral representation of the partition function and 
then use the Villain model for the three cosine terms in 
Eq. 1 [26]. Assuming that vortex rings in the 3D reser- 
voir are unimportant (which will always be valid if the 
reservoir is sufficiently dilute), we integrate out the in- 
ternal degrees of freedom of the reservoir to arrive at a 
partition function of the form 

Zvm = ^ cxp{-^^2^eU\n XT \ 2 + -^\m XT \ 

"2^8 E ^y|^m(q)+^ n n(q)| 2 }. (2) 

Here (3 = (fcT) -1 — * 00, L x is the number of sites in the 
ID chain, and e is the time slice interval. Physically, the 
integer variables m XT may be understood as bond cur- 
rents and the fluctuations in the site occupation 
number. Eq. 2 is a form of the Bose-Hubbard model. The 
coefficient l/h(q) contains information about the reser- 
voir degrees of freedom, in particular the gapless collec- 
tive mode it supports due to its own superfluidity, 
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where cr = ^/sHrUb, A is the momentum cut-off, and 
7 = (l/47rei fl ) [27]. 

The partition function Zvm = ^ m „exp(-H™) 
may be reexpressed in terms of another pair of integer 
fields <P(x,t) and A(x,t) with m(x,r) = — <9 T 0(x,r) and 
n(x, t) = d x (j)(x, t) + A(x, t), so that 
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Here, K = et and c = veHJt. Configurations for which 
V(/> 7^ and A = contain closed loops which may be 
understood as worldlines of particle-hole pairs that sep- 
arate and recombine. The A field, which should be re- 
garded as residing on the time interval links, can cancel 
the gradient energy (dxip) 2 on part of a closed loop con- 
figuration to form individual particle or hole trajectories; 
the endpoints (occuring where d T A(x,T) 7^ 0) represent 
tunneling events between the 3D and ID systems. Tun- 
neling events may be shown [20] to be dual to the vortices. 
Moreover, the model defined by Eq. 4 may be generalized 
to include a core energy E c for the vortices. 

As we discuss below, this generalized Hamiltonian sup- 
ports three phases. When tunneling events proliferate 
through the system, the ID and 3D systems exchange 
particles freely (equivalently, vortices of the 1+1 dimen- 
sional system are linearly confined), and their phases be- 
come locked together to form a single superfluid. As 
fluctuations in <f> and A are decreased (by decreasing E c 
and/or h, or decreasing K), tunneling events bind into 
pairs which conserve the overall particle number in the 
ID system, although closed particle-hole worldlines re- 
main proliferated. This state may be understood as one 
in which Josephson vortices between the 3D and ID sys- 
tem proliferate due to quantum fluctuations, effectively 
decoupling the two systems. For still smaller fluctuations, 
particle-hole worldlinc loops of arbitrarily large size be- 
come irrelevant, indicating that the fluctuations in the 
particle number on each site have been suppressed, and 
the system is a Mott insulator. 

RG analysis: A method for constructing a momen- 
tum shell RG for Hamiltonians such as Eq. 4 was devel- 
oped in Ref. [20]. We replace the integer fields <j> and 
A with the continuous fields (p(x, r) and a(x, r) in Eq. 4 
and add terms of the form — y f dxdr cos{2it<p(x,t)] and 
— y a J dxdr cos[27ra(x, r)] so that the resulting effective 
Hamiltonian has the same symmetries as the original 
one [21]. In this replacement, y — cxp(— E c ) is the usual 
vortex fugacity, and choosing y 2 ^ J d 2 ql/h(q) approx- 
imately reproduces the action associated with a world- 
line endpoint (i.e., a tunneling event). We then inte- 
grate out short wavelength degrees of freedom [A/6 < 
\Qx\, \w n \/c < A with b = exp(Z)] to lowest order in 
y and y a , and rescale lengths, times, and fields ac- 
cording to x = bx',T = br' , if' (x' , r') = <f(x,r), and 



a(x,r) = a'(x' ,r')/b. This choice preserves the terms 
that are lowest order in gradients in the quadratic part of 
the Hamiltonian. Because the a field shrinks upon rescal- 
ing, it is natural to expand cos[27ra(x, r)] in its argument, 
producing a quadratic term of the form ^p\a(x, t)| 2 that 
contributes to the fixed point. The higher order vertices 
generated by this expansion contribute to the renormal- 
ization of p but are themselves irrelevant. Note the initial 
value of p is A-K 2 y a . 

The fixed points that emerge from this procedure have 
the form 
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The last term in is very important: since tunneling 
events are specified by d T a ^ 0, when p ^ they are 
bound into equal and opposite pairs. Unbinding occurs 
if p{tj scales to zero. To lowest order in y a , its scaling 
relation is 
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For small p one can easily show that the term 
—y J dxdr cos[2ir ip (x, t)] is strongly irrelevant. The re- 
sulting scaling flows are shown in Fig. 1. As may be 
seen, in general p scales to a point along a fixed line; if 
its initial value is small enough (as occurs for small 7 
and large J), that point is at p — 0. A remarkable fea- 
ture of Eq. 6 for such flows is that p(£) = at a finite 
value of I — £*, because of the singularity as p — > 0. If 
the initial value of p is increased, I* eventually diverges, 
defining a transition point above which p scales to a non- 
zero value. This represents a confinement transition for 
tunneling events [20,28]. 

The dual representation of this system is obtained by 
performing a Poisson resummation on A and <j> in Eq. 4. 
The resulting Hamiltonian is very similar in form to 
Eq. 4 [29] , and may be understood as a model with vortex 
degrees of freedom rather than tunneling events. Because 
of this, there is a second transition, in which vortices go 
from a bound to an unbound state, in the same continu- 
ous fashion as seen above for tunneling events. In terms 
of f and a degrees of freedom, this corresponds to a phase 
in which arbitrarily large particle-hole loops are negligi- 
ble in the partition function. Pairs of tunneling events 
are then linearly confined [20]. 

The phases described above have physically different 
characteristics. This can be seen most directly in the col- 
lective mode spectra, which in principle may be measured 
in light scattering experiments [30,31]. We find these 
by examining the density-density correlation function 
(n(— q x , — ui)n(q x , u))) near the fixed point Hamiltonians 
representing the various phases. Generally, this contains 
a broad response for lo > cnq x due to modes in the reser- 
voir. Beyond this, in the superfuid phase [Fig. 2(a)], 
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we find [30] a sharp resonance (collective mode) at uoi — 
c R q x -8uj, where Suj = c R Aexp{-(c 2 /-fK ) / {ql[c 2 R - c 2 ])} . 
This rapidly approaches the reservoir continuum as q x — > 
0, indicating that number fluctuations in the ID system 
will strongly mix with those of the reservoir in the long 
wavelength limit. This is consistent with our interpreta- 
tion of this phase as a single 3D superfluid. (Note that 
although the term containing the logarithmic singularity 
is irrelevant in the RG sense, it nevertheless has the phys- 
ical effect of "pulling" the collective mode very close to 
the edge of the reservoir modes at small enough wavevec- 
tors.) In the intermediate phase we find two collective 
modes, the superfluid mode of the reservoir at w\ and 
another linear mode at 102 = cq x pK / (pK + c 2 ) [Fig. 
2(b)]. For an appropriate geometry we will show this 
leads to metallic behavior. In the Mott insulator phase, 
the reservoir becomes decoupled from the ID chain in 
the long- wavelength limit, and supports a gapped mode 
at co 3 = ^c 2 q 2 x + 2tt 2 K/E c [Fig. 2(c)]. 

Conductance: The different characters of the phases 
may also be seen in the conductance [32] of the system. 
To measure this, one has to attach a particle source and 
lead as shown in Fig. 3. The conductance quantifies the 
current injected by this source with chemical potential 
p s > 0, draining into the 3D reservoir which is held at 
zero chemical potential. (Note that a link has been re- 
moved to ensure that the current flows in the ID chain 
before tunneling into the reservoir.) This system could 
be fabricated in a Josephson-junction array [7]. 

The source and lead introduce two extra degrees of 
freedom \(f> L (r), N l {t)\ to our model. Wc treat the par- 
ticle source as an ideal reservoir with Hamiltonian Hl = 
PsNl; current conservation at the site (x = 0) where 
current is injected specifies Nl in terms of the other 
variables through a constraint in the partition function, 
U T {d x m{x = 0,r) + d T n(x = 0,r) - dN L (r)/dT = 0}. 
The conductance is given by G{uo) = ui(Nl(—oj)Nl(co)). 

The analysis proceeds in a fashion similar to what is 
described above for the uniform chain. However, because 
of the lead there is an additional term in our effective 
Hamiltonian of the form — Sy J dr cos[27ra(a; = 0, r)], re- 
flecting the fact that the rate of tunneling events at the 
x = site is different than at other sites. The scaling 
relation for 5y(£) takes the form [30] 



necessarily grow, and we need to look for a new fixed 
point. To do this, we integrate the RG flows to a scale £q 
for which the irrelevant operators may be ignored. We 
then integrate out if and a(x ^ 0, r) and recollect some 
irrelevant terms to restore the ^pa(x = 0,t) 2 term to 
its cosine form cos[27ra(£o)a(:c = 0, r)], and arrive at an 
effective Hamiltonian for the lead site in the chain, 
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We then modify the RG so that a'{x = 0,t') = a(x = 
0, r) to preserve the form of the quadratic term in "H^jp 
the scaling relation obeyed by yo is then 
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It is important to recognize that a(£o) does not shrink 
as it did in Eq. 7. Thus since p is small, yo is irrelevant. 
Our fixed point is then 7~Leff with 2/0 = 0. In this case 
the conductance is finite as uj — > 0, with G(u>) cx 1/^/jo. 
This metallic behavior is surprising in light of the linear 
mode supported by this phase. It is a result of the very 
limited phase space available for fluctuations in one di- 
mension [32]. For this reason the transport properties 
of this intermediate phase is distinct from that of the 3D 
superfluid phase. 

The conductance in the deconfincd vortex state can 
be computed straightforwardly using the dual represen- 
tation of the model. The result unsurprisingly is an in- 
sulating response G(uu) <~ u. (Details will be presented 
elsewhere [30].) 

In summary, we have shown that bosons in a one- 
dimensional optical lattice that exchanges particles with 
a bulk superfluid supports three distinct states, with dif- 
ferent collective mode spectra and conductances. 

The authors thank K. MacAdam for helpful discus- 
sions. Support was provided by the NSF under Grant 
Nos. PHY9907949, DMR0108451, and DMR0454699. 
HAF thanks the KITP for its hospitality. 
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for small y, y a , and Sy, with a{£) = cxp(—£). It is impor- 
tant to recognize that in the 3D superfluid state, p{£) — ► 
for finite £ = £*, so that 5y(£) will be driven to zero at 
1 = 1*. Thus the non- uniformity of tunneling events 
along the chain is irrelevant in this state. In this situ- 
ation, we can compute the conductance using the fixed 
point Hamiltonian (Eq. 5) for p = to find a true super- 
conducting response, G(ui) cx —i/(u> + iS). 

By contrast, if p remains finite as £ — > 00, Sy will 
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FIG. 1. Schematic diagram of RG flows for the scaling re- 
lations of p and y from initial microscopic values to a fixed 
line, p — represents a deconfined phase for tunneling events. 
For p > they are bound in pairs such that the net particle 
flow from the reservoir vanishes. Heavy line separates flows 
for the two kinds of states. 
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FfG. 2. Collective mode spectra obtained from fixed point 
Hamiltonians. In the 3D superfluid phase there is a single 
sharp mode (a) . The ID metallic phase supports two gapless 
modes (a) and (b). In the Mott insulator phase, there is a only 
single gapped mode (c). Note all three phases also contain a 
continuum of modes due to the reservoir (shaded region). 
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FIG. 3. Josephson junction realization of the system with 
source and drain to measure conductances of the phases. 
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